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ABSTRACT 

Aiming to learn about the chemistry of the dense PDR around the ultracom- 
pact (UC) H II region in Mon R2, we have observed a series of mm-wavelength 
transitions of C3H2 and C2H. In addition, we have traced the distribution of other 
molecules, such as H^^CO^, SiO, HCO, and HC3N. These data, together with 
the reactive ions recently detected, have been considered to determine the phys- 
ical conditions and to model the PDR chemistry. We then identified two kind 
of molecules. The first group, formed by the reactive ions (C0+, HOC+) and 
small hydrocarbons (C2H, C3H2), traces the surface layers of the PDR and is pre- 
sumably exposed to a high UV field (hence we called it as ^^high UV\ or HUV). 
HUV species is expected to dominate for visual absorptions 2 < Ay < 5 mag. 
A second group (less exposed to the UV field, and hence called ''low UV\ or 
LUV) includes HCO and SiO, and is mainly present at the edges of the PDR 
{Ay > 5 mag). While the abundances of the HUV molecules can be explained 
by gas phase models, this is not the case for the studied LUV ones. Although 
some efficient gas-phase reactions might be lacking, grain chemistry sounds like 
a probable mechanism able to explain the observed enhancement of HCO and 
SiO. Within this scenario, the interaction of UV photons with grains produces an 
important effect on the molecular gas chemistry and constitutes the first evidence 
of an ionization front created by the UC H II region carving its host molecular 
cloud. The physical conditions and kinematics of the gas layer which surrounds 
the UC H II region were derived from the HUV molecules. Molecular hydrogen 
densities > 4 10^ cm~^ are required to reproduce the observations. Such high den- 
sities suggest that the H II region could be pressure- confined by the surrounding 
high density molecular gas. 
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Introduction 



Mon R2 is a nearby {D — 830 pc; Herbst & Racine 1976) complex star forming region. It 
hosts an ultracompact (UC) Hll region near its center, powered by the infrared source Mon 
R2 IRS 1. The molecular content of this region has been the subject of several observational 

studies in the last decade. So far, it is known a huge CO bipolar outflow (Meyers-Rice & 
Lada 1991), ~ 15' long (= 3.6 pc), which is probably a relic of the formation of the BOV star 
associated to IRSl (Massi, FeUi, & Simon 1985; Henning, Chini, & Pfau 1992). 

Mon R2 is an excellent laboratory to learn about the physical conditions and kinematics 
of an extreme PDR (Gq = 5 10^ in units of Habing field, Rizzo et al. 2003). Continuum 
observations at 6cm, reported by Wood & Churchwell (1989), show that the UC Hll region 
is highly asymmetric, has a cometary shape and reaches its maximum toward its exciting 
star, Mon R2 IRSl. The CS 7^6 emission map from Choi et al. (2000) shows that the 
UC Hll region is located inside a cavity and bound by a dense molecular ridge (see Fig. 1). 
This picture is consistent with infrared images of the region, which show an almost spherical 
distribution of hot dust surrounded by the molecular cloud. 

The physical and chemical conditions of the PDRs associated to UC Hll regions remain 
not enough known, mainly due to important observational problems: (a) These PDRs are 
too thin to be resolved by single-dish telescopes; (b) They are usually located in very com- 
plex star-forming regions, where confusion with the parent molecular cloud, hot cores, and 
outflows may be considerable. The study requires specific PDR tracers which allow us to 
separate their emission from other sources. 

Reactive ions and small hydrocarbons have been proposed as abundant species, which 
may trace the hot ionized/molecular gas interface (Sternberg & Dalgarno 1995). The de- 
tection of the reactive ions CO"*" and HOC+ is almost unambiguously associated to regions 
with a high ionizing flux, either PDRs or XDRs (Fuente et al. 2003; Rizzo et al. 2003; Usero 
et al. 2004; Savage & Ziurys 2004). In the case of Mon R2, Rizzo et al. (2003) have recently 
reported the detection of the reactive ions CO^ and HOC+ towards the peak of the UC Hll 
region. In addition to chemical arguments, both the morphology and velocity profile of the 
lines suggested a link between the Hll region and these species. Both ions were detected 
towards the peak of the UC Hll regions and remain undetected in the molecular cloud. Fur- 
thermore, the velocity profile of these lines is different from those of other dense gas tracers. 
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The reactive ions have an intense emission peak at 10.5 km which is not detected in the 
other molecules observed. Therefore, the emission of this velocity component may arise in 
narrow layers of molecular gas surrounding the Hll region, where the chemistry is dominated 
by UV photons. 

Small hydrocarbons have enhanced abundances in PDRs (Fuente et al. 2003, 2005; 
Teyssier ct al. 2004; Pety et al. 2005). In particular, Teyssicr et al. (2004) have detected 
C-C3H2 and C4H towards the Horsehead nebula with an abundance one order of magnitude 
larger than that predicted by gas-phase FDR models. This is clearly seen by comparing the 
[c-C3H2]/[HC3N] ratio in PDRs and dark clouds. While both species have similar abundances 
in dark clouds, the [c-C3H2]/[HC3N] ratio is above 10 in PDRs. Since both molecules arc 
easily destroyed by photodissociation, this suggests the existence of an additional C-C3H2 
formation mechanism in PDRs. 

Enhanced abundances of some other compounds have also been related to the presence 
of intense UV fields into molecular clouds. Schenewerk et al. (1988) measured a value of 
[HCO]/[H^3Qo+]=9.7 in the H II region NGC2024. However, the same authors report 
values significantly lower than 1 -i.e., an order magnitude below those found in NGC2024- 
in Galactic clouds without developed H II regions, or having no indication of star formation. 
More recently, Schilke et al. (2001) have searched for HCO in a reduced sample of prototypical 
PDRs; the estimated [HCO]/[H^^CO+] abundance ratios range from 30 (in the Orion bar) 
to 3 (in NGC7023). The largest HCO abundances arc found in the Orion bar, the paradigm 
of interaction between an H II region (M42) and its parent molecular cloud. Garcia-Burillo 
et al. (2002) obtained a high angular resolution image showing widespread HCO emission in 
the external galaxy M82. The enhanced HCO abundance ([HC0]/[H^3C0+] ~ 3.6) measured 
across the whole M82 disk was also interpreted in terms of a giant PDR of 650 pc size. 

SiO is known to be a privileged tracer of large-scale shocks in the interstellar medium 
(Martin-Pintado ct al. 1997; Garcia-Burillo et al. 2000, 2001). Its fractional abundance is 
enhanced by more than 3 orders of magnitude in shocked regions relative to its abundance 
in quiescent clouds where it takes values < 10~^^. Schilke et al. (2001) observed SiO towards 
a sample of PDRs and obtain fractional abundances of ~10~^^ in these regions, i.e., a factor 
of 10 larger than that in dark clouds. 

In this paper, we present observations of a selected set of molecular species (C2H, c- 
C3H2, HC3N, H^=^CO+, HCO, and SiO) which are thought to be good probes of the different 
layers of the molecular gas in PDRs. In this way, we have performed a study of the physical 
conditions, the chemistry, and the kinematics of the molecular gas surrounding the UC Hll 
region. 
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2. Observations 

The observations were carried out using the IRAM 30m radio telescope in Pico de Veleta 
(Spain) during July 2002 and August 2003. The multi-receiver capability of the telescope 
was used for these observations. The observed rotational transitions, frequencies, beam sizes, 
and main-beam efficiencies are shown in Table 1. The backends used for the observations 
were an autocorrelator split in several parts, a 256 x 100 kHz and two 512 x 1 MHz filter- 
banks. All the hnes but H^^CO+ 1— >-0 and C3H2 5(1,4)^4(2,3) have been measured with two 
spectral resolutions, namely the lower spectral resolution provided by the 1 MHz filter-bank 
and the higher one, 78 kHz-100 kHz, provided by the autocorrelator and the 100 kHz filter 
bank. The H^"^CO+ 1^0 and C3H2 5(1,4)^4(2,3) lines have only been measured with the 1 
MHz filter-bank, providing the lower spectral resolution. The C3H2 6(1,6) ^5(0,5) [ortho] 
and 6(0,6)^5(1,5) [para] are blended. 

The observational strategy was as follows. We have mapped the UC H II region Mon R2 
and its surroundings in millimeter transitions of H^'^CO^, HC3N, HCO, and SiO. The C2H 
and C3H2 transitions were observed along a strip transversal to the UC H II region, at an 
inclination angle (east to north) of 135°, which is approximately the symmetry axis of Mon 
R2. The mapped area (72"x72") and strip positions are sketched in Fig. 1, overlaid on the CS 
7— >6 map of Choi et al. (2000) and the 6 cm-continuum map of Wood & Churchwell (1989). 
The relative offsets are referred to Mon R2 IRSl, (RA, Dec.)2ooo=(06''07™46!2, -06°23'08'.'3). 
Longer integration time was devoted to observe two particular positions: the maximum of 
the continuum emission, at (0", 0"), and the maximum of the molecular envelope, at (+10", 
-10"). These two positions were also observed in the reactive ions CO"*" and HOC"'", and 
already published in a recent paper (Rizzo et al. 2003). 

3. Results 

3.1. C3H2 and C2H along the strip 

The small hydrocarbons C2H and C3H2 have been observed along a strip crossing the 
ionization bar in Mon R2. After the hyperfine-structure analysis -which compares the rel- 
ative intensities of each component-, we have found that the C2H triplet emission at 87.3 
GHz is optically thin at all positions. Based on this result, we will deal in the following just 
with the most intense hyperfine component (F = 2 — 1), referring to this as "the C2H 1— >0 
line". Some of the C2H and C3H2 observed spectra are shown in Fig. 2. The Gaussian fits 
to a single component are also shown in Table 2 for all the lines. The C2H 1^0 line shows 
emission at high velocities (Vlsr < 8 kms~^ and Vlsr > 12 kms~^) that is not detected (or 
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only marginally detected) in the C3H2 lines. These velocity ranges coincide with those of 
the red-shifted and blue-shifted emission in the large-scale, now inactive, outflow associated 
with IRS 1 (Meyers- Rice & Lada 1991; Giannakopoulou et al. 1997). We have not detected 
emission at high velocity in any other of the PDR tracers studied in this paper and Rizzo et 
al. (2003). Thus, wc consider that the high velocity C2H emission is not related to the UC 
H II region but to the low density gas that share the kinematics of the large-scale outflow. 
Since we are mainly interested in the study of the PDR surrounding the UC H II region, we 
have only used the C2H emission in the range 8-12 km s~^ in our calculations. 

Within the considered velocity range, the velocities and spatial distribution of the C2H 

1— >0 and C3H2 2^1 lines are remarkable similar, which suggest that both lines trace essen- 
tially the same region. However, the emission of the C3H2 2— >-l and the C3H2 6— >5 hues 
have different kinematics and spatial distribution, which became eloquent when looking at 
the central velocities in Table 2 and the shaded areas in Fig. 2. The C3H2 2— >^1 hue is 
narrower than the 6— >5 line along the whole strip. Furthermore, the central velocity of the 

2— >1 line is lower than that of the 6^5 line towards the UC H II region. When looking at 
the areas, it is remarkable the different spatial distribution of the 2^1 and 6^5 lines. The 
2^1 emission (a low excitation line) is clearly lower at the central part of the strip, while 
the 6^5 emission (a high excitation line) is particularly intense there. 

This behavior is better explained by assuming the existence of at least two gas compo- 
nents characterized by different velocities. The blue velocity component (hereafter referred 
to as CI), emitting at velocities from 8 to 10 km s~^, is dominated by the emission of 
the 2— >1 line along the whole strip, but it is not particularly intense in the 6— >5 line. This 
component remains quite uniform along the strip. The red velocity component, with emis- 
sion from 10 to 12 km s~^and hereafter referred to as C2, has intense emission in the 6— >^5 
line but it is very weak in the 2— >^1 line. C2 is especially intense towards the inner H II 
region. The distribution of both components can be interpreted as CI associated with the 
foreground molecular cloud and C2 associated with the H II region, or the interface between 
it and the molecular cloud. 

This trend is illustrated in Fig. 3, where we plot the integrated intensity of the C3H2 
6— >-5 line (Fig. 3a) and the 6^5/2^1 line intensity ratio (hereafter referred to as R62) for 
each component (Fig. 3b). R62 remains quite constant and ~ 1 for CI along the strip. 
However, it has signiflcant variations in C2. The values of R62 are ~ 0.6 outside the UC H II 
region and > 4 towards the inner region. The Fig. 3b clearly separates both components, 
and shows that C2 is tracing highly excited molecular gas, linked to the UC H II region. 
This sphtting of the values of R62 is also in agreement with the spatial distribution of both 
components. 
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Of course, the different beams involved in both hnes affect the computation of R62. We 
could estimate the beam dilution by convolving the C3H2 6— emission of the C2 down to 
the angular resolution of the 2— >1 line. This one-dimensional smoothing was performed along 
the symmetry axis, which shows the largest variations in other molecular tracers. The results 
are shown by dashed lines in Fig. 3. As expected, the distribution of the 6^5 line intensity 
peaks in the southeastern border of the UC H II region. Furthermore, it is remarkable that 
R62 remains > 5 in the three central positions. In other words, the high value of there, 
even smoothed, strengthens the mutual relationship between C2 and the UC H II region. 

3.2. Maps 

Wc have mapped the region indicated in Fig. 1 in several rotational lines. Fig. 4 shows 
the resultant maps in Hi3C0+ 1^0 (Fig. 4a), SiO 2^1(Fig. 4b), HCO 1^0 (Fig. 4c), 
and HC3N 10^9 (Fig. 4d). Superimposed are the 6 cm-continuum emission from Wood & 
Churchwell (1989). All the maps were constructed by integrating the lines along the total 
velocity range of emission. We note remarkable differences in the distribution of the different 
species. The H^^CO''" and the HC3N emission clearly follow the high density gas, properly 
traced by the CS 7^6 emission (see Fig. 1) and reach a minimum toward the center of the 
UC H II region. This is consistent with the interpretation of H^^CO"'" and HC3N as being 
photodissociated within the UC H II region. 

However, the other molecules present different behaviors. The SiO and HCO emissions 
seem to be located preferably at intermediate positions between the molecular envelope and 
the H II region. Although this effect should be regarded as tentative, because the spatial 
difference between the continuum peak and the SiO distribution is about half a beam, the 
systematic difference looks evident, and clearly the SiO or HCO maps are different to those of 
H^^CO+ or HC3N. This result is in line with Schilke et al. (2001), who also found moderately 
enhanced SiO and HCO emissions close to the ionization front in the prototypical H II region 
M42. However, the two species peak at different positions. The HCO maxima appears almost 
exactly at the SiO minimum, and both emissions completely bound the inner UC H II region. 

In Fig. 5 we show all the spectra we have observed towards the (0",0") position (this 
paper and Rizzo et al. 2003). The dashed vertical hne traces the hmits of CI and C2. 
Unfortunately, some of the lines at (0",0") are rather weak, which yields a low S/N ratio. 
Even though, a quick look suggests the existence of two different behaviors in the shape of 
the observed lines. A first group of lines has the maximum intensity at velocities < 10 km 
s~^, close to the CI velocity; this is the case for the SiO and HCO lines. A second group 
peaks at velocities > 10 km s~^ and are especially intense in C2; the C3H2 6^5, CO"*" 2^1 
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and HOC+ 1— >0 belongs to this group. The strong difference between the SiO hnes at C2 
velocities should be more carefully analyzed in future observations. 

After a joint view of the maps (Fig. 4) and the spectra (Fig. 5), we see that the SiO and 
HCO emitting gas do not seem to share the same volume with the reactive ions and small 
hydrocarbons, but to arise in a more shielded layer of the PDR, which is preferably at the 
velocity of CI. As further discussed in Sects. 5 and 6, different excitation conditions and 
chemistry are associated to this shielded layer and C2. 

4. Physical conditions: LVG results for the hydrocarbons 

We have developed a LVG code to estimate the main physical parameters of the C3H2 
molecule. The Einstein A- coefficients have been taken from Chandra (2003). The coUisional 

rates are from Chandra & Kegel (2000), who provided values for kinetic temperatures (Tk) of 
30, 60, 90, and 120 K. Our LVG code treats the ortho- and para-C3H2 as different species. In 
most positions, we have detected only the 2(1,2)^1(0,1) [ortho] and the 6^5 (a blending of 
the 6(1,6)^5(0,5) [ortho] and 6(0,6)^5(1,5) [para]) lines. Then, we need to assume a value 
of the C3H2 ortho-to-para ratio (OTPR) in our calculations. In order to have an estimate 
of the OTPR, we have used the 6— >5/5— >4 hue intensity ratio (hereafter Res) towards the 
(0",0") and (+10", -10") positions. The observed value, Rgs ~ 1.4 ± 0.3, is compatible 
with Rq5 = 4/3, which is expected for the standard value of 3. Therefore, we have assumed 
OTPR = 3 in our LVG calculations hereafter. 

The LVG code was run for different sets of Tk, from 10 to 150 K. For Tk different to 
those tabulated by Chandra & Kegel (2000), we have interpolated or extrapolated between 
the two closest temperatures. A minimum value of 30 K for Tk is needed to reproduce the 
6— >5/2— s>l line ratio (hereafter referred to as R62) both in CI and C2. For Tk > 50 K, R62 
is weakly dependent on the assumed kinetic temperature and traces mainly the hydrogen 
density. In Table 3, we show the LVG results assuming Tk = 50 K, the kinetic temperature 
derived by Giannakopoulou et al. (1997) from multiline molecular observations. 

We have obtained rather uniform values of both the molecular density (n(Il2)) and the 
C3II2 column density (A^(C3ll2)) in CI across the observed strip. The derived n(}l2) is 
around 5 10^ cm^"^, in agreement with results derived from other tracers (Choi et al. 2000; 
Rizzo et al. 2003), and the iV(C3H2) are in the range 1-3 10^^ cm~^. A different behavior is 
observed in C2. The C2 hydrogen density increases from a few 10^ cm~^ outside the UC H II 
region to > 4 10^ cm"^ in the three central positions. This limit is almost insensitive to Tk 
in the range 50 to 150 K. Thus, the larger values of the R62 observed towards the H II region 
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cannot be due to a higher kinetic temperature. The Fig. 6 shows the LVG results from two 
of the runs, corresponding to values of Tk of 60 and 120 K. The shadowed areas indicate the 
range of observed values towards the (0",0") position. While A'^(C3H2) is similar in CI and 
C2, the C2 density is a factor > 5 greater than in CI. 

Recent results of Jaffe et al. (2003) and near-IR images show a shell-hke structure, 
intense in mid- and ncar-IR, of a few arcsec width, wrapping up the UC H II region. It is 
indeed possible that C2 is tracing the high density molecular gas associated to this feature. 
This would also explain the spatial distribution of A^(C3H2) in C2 which slightly increases 
towards the border of the H II region, the expected behavior of a molecular shell surrounding 
the ionized gas. 

We have estimated the C2H column density (A^(C2H)) across the studied strip using 
the LVG code and assuming n(H2) = 10® cm~^ and Tk = 50K. We have adopted the cross 
sections provided by Green & Thaddeus (1974) for the HCN molecule, and then used the lOS 
(Infinite-Order Sudden) approximation for molecular collision dynamics to get the coUisional 
coefficients between different levels. A similar procedure was followed by Truong-Bach et al. 
(1987), who also provided the Einstein A- coefficients for C2H. 

The CI and C2 components have been separated towards the (0",0") and (-1-10",— 10"). 
Most of the C2H emission seem to be associated to CI outwards these positions. The results 
are shown in Tables 3 and 4. A'^(C2H) varies from 2 to 4 10^^ cm~^, decreasing toward the 
center of the H II region, as expected in a region where the molecules are photodestroyed by 
a high incident UV flux. When compared to C3H2, the [C3H2]/[C2H] ratio slightly increases 
towards the center; an exception is the last point of the strip, which has a considerable lower 
S/N ratio (see Table 3). Comparing both components, the [C3H2]/[C2H] ratio is three times 
greater in C2 than in CI (sec Table 4). As we will refer in Sect. 6.1, the significantly larger 
[C3H2]/[C2H] ratio in C2 reveals that the chemistry of the molecular gas is heavily affected 
by the UV radiation. 

5. SiO and HC3N fractional abundances 

The H^'^CO^, SiO and HC3N column densities have been estimated using an LVG code 
and assuming Tk = 50 K and n(H2) = 10® cm~^. Under these physical conditions, the excita- 
tion temperature of the H^^CO"'" 1— >0 line is ~20K. We have also derived the HCO column 
density by assuming optically thin emission, LTE conditions and a rotation temperature of 
20 K. The uncertainty in the excitation temperature of HCO should not affect the fractional 
abundances of these molecules. The results for some selected positions are shown in Table 
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5. 

Significant changes in the fractional abundances are detected across the UC H II re- 
gion. In particular, the HC3N abundance decreases by a factor of ~3 towards the UC H II 
region (X(HC3N) ~ a few 10"^'^), when compared towards the dense molecular cloud (where 
X(HC3N) ~ 10~^). This decrease of the HC3N abundance in regions of enhanced UV field 

has also been found in other regions like the prototypical PDR of the Orion Bar (Rodriguez- 
Franco et al. 1998), or the starburst galaxy M82 (Fuente et al. 2005), and it is easily explained 
by photodissociation. 

SiO has a different behavior. The SiO abundance is low towards the center of the UC 
H II region (X(SiO) ~ 510~^^). It is also low towards the dense molecular cloud, where we 
derive X(SiO) < 7 10^^^, in agreement with the values of the SiO abundance measured in 
dark clouds. However, it is one order of magnitude higher at intermediate positions between 
the UC H II region and the molecular cloud; here, the SiO abundance is comparable to that 
found by Schilke et al. (2001) towards the PDRs associated with the Orion Bar and S140. 

Similarly to SiO, the HCO abundance is also enhanced at the border of the UC H II 
region. In particular the [HCO]/[H^^CO"'"] ratio is 10 times larger in the position (0",24") 
than in the center of the UC H II region and in the dense molecular cloud. We have measured 
a HCO abundance as high as ~ 10~^ towards the HCO peak. Such large values of the HCO 
abundance have been found at the border of the H II region NGC 2024 by Schenewerk et 
al. (1988) and in the Orion Bar by Schilke et al. (2001), suggesting that this molecule is 
specially abundant in the PDRs formed at the edges of the H II regions. Garcia-Burillo et 
al. (2002) found large values of the HCO abundance in the starburst galaxy M82, where a 
low-density ionized component is filling a substantial fraction of its volume (see also Seaquist 
et al. 1996). 

6. Photon-dominated chemistry in Mon R2 

6.1. The surface layers of the PDR: HUV species 

Our data show the existence of a high-density layer of gas associated to the UC H II 
region in Mon R2. This layer is traced by C2 and seems to have different physical conditions 
from that of CI and the molecular cloud. Furthermore, a different chemistry may also 
carry out in the layer traced by C2. In order to explore this idea, we have computed the 
column densities and fractional abundances for the two components, towards the positions 
(-|-10",-10") and (0",0"), and for the molecules studied in this paper and Rizzo et al. (2003). 
The results are shown in Table 4. Significant differences exist between the [00+] /[HCO"*"], 
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[HCO+]/[HOC+] and [CgHsJ/fCaH] ratios between CI and C2. In order to characterize 
the chemistry of the two components, we have modeled the PDR using the plane-parallel 
model developed by Le Bourlot and collaborators (Le Bourlot et al. 1993) with an UV field 
Go = 510^ in units of the Habing field and a total density {i.e. atomic plus molecular) of 
2 10^ cm~^ (Rizzo et al. 2003, and this paper). The main results are sketched in Fig. 7. 

The molecules CO^, HOC"'', C2H, and C3H2 arc known to be good tracers of the surface 
layers {Ay < 5 mag) of the PDRs (see e.g. Sternberg & Dalgarno 1995). These environments 
are highly exposed to an intense UV field, and hence we will refer to these molecules as "high 
UV^ species (HUV for briefing). In particular, the reactive ion CO^ is only expected to reach 
significant abundances at a visual extinction. Ay < 2 mag (Fuente et al. 2003; Rizzo et al. 
2003). In Fig. 7a we show the C0+ and HCO+ fractional abundances as a function of the 
visual extinction for our model. As expected, the abundance of CO"*" is only significant 
at Ay < 2 mag, and just at this low visual extinction, the predicted [CO"'"]/[HCO"'"] ratio 
agrees with that measured in C2 towards the H 11 region (light gray region in Fig. 7a). 
The [CO"'"]/[HCO+] ratio is more than one order of magnitude greater in C2 than in CI 
(Table 4), which indicates that CI arise in a more shielded layer of the molecular gas having 
Ay > 4 mag (dark gray region in Fig. 7a) . 

The [HCO"'"]/[HOC"'"] ratio is a factor of 2 lower in C2 than in CI. Low values of the 
[HCO^]/[HOC+] ratio are associated to highly ionized regions (Usero et al. 2004; Savage & 
Ziurys 2004; Fuente et al. 2005). Unfortunately, the reactive ion HOC+ is not included in Le 
Bourlot's model, but we can estimate its abundance in an indirect way. The [HCO^]/[HOC"'"] 
ratio is strongly dependent on the electron abundance. In order to have [HCO"^]/[HOC"'"] 
~ 450, it is required an electron abundance, X(e"') > 210""^ (Usero et al. 2004). Using 
our model, we conclude that these high electron abundances are only reached at Ay < 
4 mag. Thus, the [CO-^J/HCO"^] and the [HCO^-J/fHOC"^] ratios support the interpretation 
of C2 as arising from a layer of molecular gas located at a visual extinction < 4 mag from 
the advancing ionization front. A new PDR model including reactions of HOC^ is under 
development, and it will directly predict the [HCO^]/[HOC"'"] ratio, as a function of global 
physical conditions (UV field, density, extinction) in a near future. 

The [C3H2]/[C2H] ratio also changes between CI and C2. In fact, this is a factor of 
2-3 larger in C2 than in CI. Our chemical model predicts that the [C3H2]/[C2H] ratio is 
maximum at a visual extinction between 2 and 4 mag (see Fig. 7). These values of the visual 
extinction are in agreement with those derived from the measured [IICO"'"]/[HOC"'"] ratio. 
Thus, it seems reasonable that the emission of C2 arises in the same layer of dense gas 
than HOC"*" {Ay < 4 mag), while CO"*" arises in even a more external layer {Ay < 2 mag). 
However, our model fails to predict values of [C3H2]/[C2H] as large as those measured in C2 



- 11 - 



even at this low visual extinction. The maximum value of the [C3H2]/[C2H] ratio predicted 

by the model is ~ 0.01, three times lower than those measured in C2. (see Fig. 7b and 
Table 4). A failure of PDR models to account for the observed C3H2 abundance has already 
been commented in previous works (Fuentc ct al. 2003; Teyssier ct al. 2004; Pety et al. 
2005). The existence of an additional formation mechanism for C3H2 linked to the PAHs 
photodestruction has been proposed as an alternative to account the large C3H2 fractional 
abundances in PDRs. 

Summarizing, we can infer from the comparison of our observations of HUV molecules 
and model calculations that C2 is tracing the dense molecular gas located at a visual extinc- 
tion < 4 mag from the ionization front, while CI is mainly tracing the more shielded part 
{Ay > 5 mag) of the PDR. Consequently, the UC H II region is surrounding by a thin layer 
(< 4 10^^ cm~^) of dense (> 4 10^ cm~^) molecular gas, whose chemistry is heavily affected 
by the UV stellar radiation and kinetically well represented by C2. 

6.2. The back layers of the PDR: LUV species 

In contrast to the HUV species, the column densities of HCO and SiO decrease towards 
the UC H II region, as expected in this harsh environment. Furthermore, their abundances 
are larger in CI than in C2 (see [HC0]/[Hi3C0+] and [SiO]/[Hi3CO+] in Table 4), and show 
an enhancement at its edges (see Table 5). 

Clearly, both HCO and SiO do not share a behavior similar to the HUV molecules; 
moreover, they are not particularly abundant in the cold envelope. Both molecules seem 
to be mainly present in the shielded parts of the PDR (^v > 5 mag) , and would belong 
to a family of species which can survive in those environments, less exposed to the UV 
field than the HUV molecules. We hereafter call such molecules as "/ow; C/V" (or LUV for 
briefing) species. The abundance of our LUV molecules are unusually high and could not be 
explained by gas-phase models (see, for example, Leung ct al. 1984; Sternberg & Dalgarno 
1995). Equally puzzling is the fact that both molecules peak at different positions. 

As we have already stated before, both SiO and HCO have -near their maxima- abun- 
dances comparable to other PDRs, such as S140 and Orion Bar (Schilke et al. 2001). More- 
over, the [HCO]/[H^^CO^] abundance ratio near the HCO peak is in the range 10-15, higher 
than in NGC2024 (Schenewerk et al. 1988) and M82 (Garcia-Burillo et al. 2002). This en- 
hancement can be due to a few reasons. Firstly, some relevant chemical reactions leading 
to formation of SiO and HCO might not be considered in gas-phase models. It was already 
analyzed by several authors (see Walmsley et al. 1999; Schilke et al. 2001) and sounds rather 
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unlikely. Even though, we should not disregard the hypothesis, mainly due to a possible lack 
of complete chemical networks and rate coefficients. 

So far, the most widely accepted alternative explanation arises by considering the grain 
chemistry. As many authors pointed out (see, for example, the review by Caselli 2005), the 
dust plays a key role in driving the chemistry and the energetics in PDRs. Walmsley et al. 

(1999) have modeled the Si chemistry in PDRs, and found that most of Si is in solid form 
at extinctions greater than 3 mag. This issue was observationally confirmed by Fuente et 
al. (2000) in NGC 7023, where a special distribution and enhancement of SiO was detected. 
Schilke et al. (2001) have included photodesorption of Si in their FDR model. SiO would later 
be produced in gaseous phase by reactions with OH and O2. By doing so, those authors could 
account for the observed SiO abundance in Orion Bar and S140. Laboratory measurements 
of the rate coefficient of the reaction Si + O2 — > SiO + O (Le Picard et al. 2001), followed 
by detailed calculations of Si chemistry, have confirmed the results by Schilke et al. (2001). 

In the case of HCO, the photodesorption of H2CO -followed by its photodissociation in 
gaseous phase- have been claimed as a mechanism capable to counterbalance the photodis- 
sociation of HCO against the UV photons (Schilke et al. 2001). The authors could get an 
enhancement of HCO (at Ay > 5 mag) but not to the observed levels. This may be due to 
the need of more reliable rates for the processes involved, but also to the reduced network 
used in the gas-phase formation of this low-abundance species. We have measured similar 
abundances of HCO and hence are in a similar situation; the UV field intensity, greater in 
Mon R2 than in Orion Bar (Rizzo et al. 2003), does not seem to affect the HCO abundance. 

However, photodesorption may not efficiently work in some other environments. Ruffie 
& Herbst (2001) have modeled both gas- and dust-chemistry in quiescent cores, and shown 
that the abundance of some species with lower binding energies than Si do not significantly 
change when including photodesorption. The modeled quiescent cores, however, have a Go/n 
ratio at least two orders of magnitude lower than in Mon R2, because the densities are two 
orders of magnitude lower, and the incident UV field takes interstellar values. Even though, 
their results should be taken into account, and a complete model like the Ruffie & Herbst's 
one, under physical conditions more typical of PDRs, will surely improve our knowledge 
about the particular chemistry of LUV species and PDRs in general. Anyway, if SiO and 
HCO are formed after photodesorption of Si and H2CO, respectively, we would expect other 
atoms and molecules do so, and would be detectable. At this respect, an interesting test 
may be to carry out an observational campaign to search for and study other LUV species 
throughout the whole field. A good candidate to belong to the LUV group in NH3, a widely 
recognized species linked to grain chemistry, whose abundance close to the borders of H II 
regions is enhanced (Larsson et al. 2003). 
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The different spatial distributions of our LUV molecules is another interesting issue. 
Both molecules almost complete a "ring" between the UC H II region and the molecular 
envelope. However, they peak at different locations throughout the ring: while SiO arises 
in the densest part of the ionization front, HCO peaks in the tail of the cometary shape. 
Different local conditions should exist in order to explain the variations in relative abun- 
dances of SiO and HCO. A rather simple explanation is that HCO and SiO have a different 
desorption yield; in this case one would expect a layered structure in the morphology of these 
molecules. Unfortunately, the low angular resolution of the HCO and SiO maps prevent us 
from concluding about that. A second possibility is that a low- velocity shock associated 
with the ionization front contributes to enhance the HCO and/or SiO abundances. In the 
case of hot cores, this process was already analyzed by Viti et al. (2001). The low velocity 
shocks contribute to release Si and/or H2CO from the grain surfaces, and the gas phase 
chemistry is modified by the advancing shock front. The radical HCO is rapidly destroyed 
by O, but in shocked regions the O abundance is reduced because it is converted to H2O. 
Furthermore, SiO is efficiently formed by the reaction Si + OH SiO + H, being OH the 
photodissociation product of H2O. Thus, the HCO and SiO abundances are very sensitive 
to the O/OH/H2O relative fractional abundances and can change locally because of low ve- 
locity shocks and small inhomogeneities in the incident UV field. Finally, a variation in the 
grain composition may account for the differences between the HCO-peak to the SiO-peak. 

Observations of dust emission may help to solve this interesting puzzle in Mon R2. Krae- 
mer et al. (2001) have studied this region in the mid-IR, from 8.2 to 20.6 iim. Besides the IR 
sources, continuum dust emission is dominated by a ring-like structure which roughly bounds 
the ionized region. While the dust temperature is rather uniform in the area (100-115 K), 
the opacity is greater in the ring. Kracmcr et al. (2001) also studied the silicate absorption 
feature at 10 /xm, and found a patchy distribution, not necessarily correlated to the point 
sources, the temperature, or the opacity. So far, dust properties may significantly change 
in areas as small as a few arcsec. Although this result encourages one of our alternatives 
(variation in grain composition), the distribution of HCO and SiO is still an open subject 
because the other proposed explanations (different desorption yield and low- velocity shocks) 
may be working in this case. 

To summarize, our SiO and HCO observations show that these compounds are enhanced 
in the outskirts of the UC H II regions (^v > 5 mag) . The enhancement in the abundance of 
these species may be related to the interaction of the UV photons with grains, instead of the 
direct effect of the UV photons on the molecular gas chemistry. The interaction of the UV 
photons with grains may modify the chemistry of large column densities of the molecular 
gas, and may show an evidence of the advancing PDR into the molecular cloud. 



-14- 



7. On the kinematics of the dense gas layer 

In the UC H II regions, the high densities of the ionized gas lead to enormous internal 
pressures {nT ~ 10^ cm~^ K). Since the sources are very small, their expansion timescales 
should be as short as a few hundred years (Dreher & Welch 1981). Surprisingly, the number 
of UC H II regions is greater than the expected from the massive star formation rate (Wood 
& Churchwell 1989), which is known as the "lifetime problem" of the UC H II regions. A 
deep knowledge of the morphology and physical conditions of the dense PDR around the 
UC H II region in Mon R2, would shed some light on this problem. We can explore one of 
the proposed explanations, in which the extremely dense neutral gas confine the Hll region 
by pressure (de Free et al. 1995; Akeson & Carlstrom 1996); our data show the existence of 
a gas layer with n(H2)> 4 10^ cm^^. Hence, if we assume a kinetic temperature of 100 K, a 
hydrogen density n{B.2) ~ 10^ cm^'^ would be enough to pressure-confine the UC H II region. 

The velocity range associated to C2 is undoubtedly different from that linked to the 
parent molecular cloud, and hence it should have a particular kinematics. By means of the 
PDR chemistry, we have a good chance of further studying the kinematics of the molecular 
gas confining the UC H II region, because some chemical species may avoid the confusion of 
the foreground molecular cloud. In Fig. 8, we show the position- velocity (P-V) diagram of 
the C3H2 2— >-l (Fig. 8a) and 6— >^5 lines (Fig. 8b) across the observed strip. As we previously 
discussed, most of the C3H2 2— >^1 line traces CI, whereas the 6— >-5 emission is dominated 
by C2. The P-V diagram of the C3H2 6^5 line is preferably located toward the H II region 
(center of the map), rather than extended along the whole strip. The curves superimposed 
in Fig. 8 were symmetrically traced around CI, and roughly follow the C3H2 6^5 maxima. 
As the extreme velocities of the curves are ^ 7 and 12 km s^^, this pattern is compatible 
with a picture where the PDR is expanding at a velocity between 2 and 3 km s~^. It is worth 
mentioning that such a low expansion velocity was recently observed in an UC H II region 
in W48, by observations of carbon recombination fines (Roshi et al. 2005). 



8. Conclusions 

We have carried out a molecular survey towards the UC H II region Mon R2 to have 
a deeper insight into the chemistry associated to this kind of objects. Our results show 
the existence of two groups of molecules. The first group (HUV) is formed by molecules 
that present large abundances in the surface layers of PDRs {Ay < 5 mag). The species 
C0+, HOC+, C2H, and C3H2 belong to this group. Gas-phase PDR models can successfully 
account for the behavior of CO"*", HOC+, and C2H. In the case of C3H2, however, PDR 
models fall short of explaining the observed abundances. New C3H2 formation mechanisms 
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linked to the photodestruction of PAHs have been proposed to account for this difference. 

HCO and SiO belong to another group of molecules (LUV). Their abundances are 
enhanced at the edges of the PDR {Ay > 5 mag) , which cannot be explained by current gas- 
phase FDR models. Photodesorption from grain mantles appears like a probable mechanism 
capable of explaining the observed abundances. A revision of gas-phase chemical networks, 
as well as a model able to include photodesorption in PDRs, may significantly improve the 
knowledge of LUV species. Other molecules, like NH3, whose chemistry is related with the 
grains also present large abundances at the edges of the H II regions (Larsson et al. 2003). 
This issue may show that the interactions of the UV radiation with the grains could produce 
important effects on the chemistry of the molecular gas, even at local visual extinction as 
deep as 10 mag. 

Finally, we have used FDR chemistry to determine the physical conditions and kine- 
matics of the layer of gas surrounding the UC H II region. Using a LVG code to fit our C3H2 
observations, we have derived hydrogen densities > 4 10^ cm~^ for this layer. The proposed 
scenario of a pressure-confined H II region requires n{B.2) ~ 10^ cm~^, perfectly compatible 
with this result. In addition, the P-V diagram of the C3H2 lines might show that while the 
gas layers at a visual extinction < 5 mag are expanding at a velocity of ~ 2-3 km s"^, the 
back layers of the PDR are moving at the velocities of the foreground molecular cloud. This 
kinematics supports the interpretation that the UC H II region could be confined by a dense 
(n(H2) > 10^ cm~^) and thin (< 4 10^^ cm~^) layer of neutral gas. 

On the other hand, the search for other LUV molecules looks as an exciting observational 

task which may lead to determine the chemical processes carrying out in the more shielded 
part of PDRs. Moreover, further observations and multi-transitional studies of the HUV 
molecules are needed to determine accurately the physical conditions and dynamics of this 
neutral layer. 

We are grateful to the technical staff in Pico de Velcta for their professional support 
during the observations. We specially wish to thank the anonymous referee, who greatly 
improved the paper by sharing with us his/her viewpoints about the PDR chemistry. This 
paper has been partially funded by the Spanish MCyT under projects DGES/AYA2000-927, 
ESP2001-4519-PE, ESP2002-01693, and AYA2003-06473. 
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Fig. 1. — Location of the positions observed in this paper. Greyscale represent the CS 7^6 
hne emission map from Choi et ah (2000). The circles across the symmetry axis of the UC 
Hll region indicate the positions observed in the C3H2 and C2H hnes. The small square, 
72" X 72" in size, traces the mapped area in HC3N, H^=^CO+, SiO, and HCO. 
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Fig. 2. — A sample of C3H2 6^5 (a), C3H2 2^1 (b), and C2H 1^0 (c) spectra along the 
strip. Tlie positions are indicated at the top part of the figure. The shadowed areas indicate 
the approximate velocity extension of the two kinematic component present in the region. 
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Fig. 3. — C3H2 observed parameters along the strip, for each of the velocity components 
defined in the text, (a) C3H2 6^5 line intensity, (b) C3H2 6— >5 to 2— >1 line ratio. Note 
the uniformity of the CI component and the separation of both components in Fig. b. The 
dashed curves are the corresponding to the C2 component, after smoothing the 6— >5 line to 
the resolution of the 2— >1 line. Angular distances are referred to Mon R2 IRS 1, increasing 
toward the southeast. 
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Fig. 4. — Emission maps around the UC Hll region Mon R2. The mapped rotational hnes 
are indicated in the top left corner of each map. Contour levels are 30% to 90% of the 

map peak, in steps of 10%. The map peaks are 2.265, 0.243, 0.244, and 3.651 Kkm s~^ for 
maps a to d, respectively. Superimposed are the 6 cm-continuum emission from Wood & 
Churchwell (1989). 
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Fig. 5. — Spectra toward the position (0", 0"). The rotational transition observed is indi- 
cated. Temperature scale is Tmb, in K. There are two types of emission depending on the 
velocity range, roughly in agreement with the velocity components defined in Sect. 3. The 
emission below 10 km s~^ is associated to the molecular envelope (component 1), while the 
emission above 10 km s~^ is linked to the ultracompact H II region (component 2). 
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Fig. 6. — LVG results at the (0",0") position, for kinetic temperatures of (a) 60 K and (b) 
120 K. For each component, it is plotted the observed values of the C3H2 6^5 line intensity 
(vertical/diagonal lines), as well as the C3H2 [6^5]/ [2—^1] line ratio (horizontal lines), as 
functions of the C3H2 column density and the H2 density. Shadowed areas take into account 
the observational errors. Although both components have comparable column densities, the 
CI density is lower than that corresponding to C2 at least by a factor of four. 
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Fig. 7. — Model predictions for the abundances of various species derived from updated Le 
Bourlot et al.'s model. The calculations have been carried out for n(H2)=2 10^ cm"^ and 
Go=5 10^ in units of the Habing field, (a) C0+ and HCO"*" fractional abundances. C0+ is 
significantly abundant at visual extinctions around 2 mag, indicated by the light gray area, 
which roughly corresponds to the values measured in C2. On the other hand, the dark gray 
area agrees with the upper limits of the [CO^/HCO^] ratio measured in CI. (b) C3H2 and 
C2H fractional abundances, as a function of the visual extinction. The [C3H2]/[C2H] ratio 
peaks between 2 and 5 mag (gray areas) . 
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Fig. 8. — Position-velocity diagram of (a) C3H2 2^1 and (b) C3H2 6^5 lines across the 
strip. Angular distances are referred to Mon R2 IRSl, increasing toward the southeast. 
Most of the C3H2 2^1 emission arises from the CI component, between 9 and 10 km s^^. 
The curves roughly follow the maxima at high velocities, which includes the C2 component. 
This pattern is compatible with an expansion of 2-3 km s~^ 
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Table 1. Observing frequencies and telescope parameters 



Molecule 


Transition 


Preq (GHz) 


beam 


Vmb 


C2H 


Ar = 1 - 0, J = 3/2 - 1/2 










F = 1 ^ 1 


87.284156 


28.5" 


0.77 




F = 2 ^ 1 


87.316925 


28.5" 


0.77 




F = 1 ^ 


87.328624 


28.5" 


0.77 



C,H 



3^2 





i^c) = 








2(1,2) - 


- 1(0,1) 


85.338893 


29" 


0.77 


3(2,2) - 


-2(1,1) 


155.518295 


17" 


0.65 


6(0,6) - 


- 5(1,5) 


217.822148 


11.5" 


0.52 


6(1,6) - 


- 5(0,5) 


217.822148 


11.5" 


0.52 


5(1,4) - 


- 4(2,3) 


217.940046 


11.5" 


0.52 



HC3N J =10 ^ 9 90.9789933 27" 0.76 

HCO Nk-k+ = loi - Ooo 

J = 3/2 - 1/2, F = 2 ^ 1 86.6708200 28.5" 0.77 

SiO J = 2 ^ 1 86.8468910 28.5" 0.77 

SiO J = 3^ 2 130.268702 19" 0.70 

J = 1 ^ 86.7543300 28.5" 0.77 



-28- 



Table 2. One-component Gaussian fits to the C3H2 and C2H spectral lines 



Offset 


Line 


Area 


Vlsr 


AVi/2 


arcsec 




K km s~^ 


km s~^ 


km s~^ 



(20, 


-20) 


C3H2 


6 - 


-^5^ 


< 0.13 










(10, 


-10) 


C3H2 


6 - 


->5^ 


1.10 (0.10) 


10.9 


(0.1) 


2.8 


(0.3) 


(0 


,0) 


C3H2 


6 - 


->5^ 


0.82 (0.11) 


10.4 


(0.2) 


3.7 


(0.6) 


(-10, 


+10) 


C3H2 


6 - 


5^ 


0.76 (0.10) 


9.8 


(0.2) 


3.6 


(0.6) 


(-20, 


+20) 


C3H2 


6 - 


5"" 


0.76 (0.15) 


10.5 


(0.3) 


3.3 


(0.8) 


(-30, 


+30) 


C3H2 


6 - 


5^ 


0.70 (0.14) 


9.5 


(0.3) 


2.7 


(0.5) 



(20, -20) 


C3H2 5(1,4) - 


-^4(2, 


3) 


< 0.07 








(10, -10) 


C3H2 5(1,4) - 


-^4(2, 


3) 


0.76 (0.09) 


10.6 


(0.3) 


3.4 (0.5) 


(0, 0) 


C3H2 5(1,4) - 


4(2, 


3) 


0.56 (0.07) 


8.9 


(0.3) 


4.2 (0.6) 


(-10, +10) 


C3H2 5(1,4) - 


4(2, 


3) 


0.47 (0.12) 


10.6 


(0.6) 


4.4 (1.2) 


(-20, +20) 


C3H2 5(1,4) - 


-^4(2, 


3) 


< 0.06 








(-30, +30) 


C3H2 5(1, 4) - 


-^4(2, 


3) 


< 0.06 








(20, -20) 


C3H2 3(2, 2) - 


-> 2(1, 


1) 


< 0.11 








(10, -10) 


C3H2 3(2, 2) - 


2(1, 


1) 


0.18 (0.04) 


11.6 


(0.1) 


0.5 (0.2) 


(0, 0) 


C3H2 3(2,2) - 


2(1, 


1) 


0.38 (0.09) 


9.8 


(0.7) 


5.6 (1.3) 


(-10, +10) 


C3H2 3(2, 2) - 


2(1, 


1) 


< 0.09 








(-20, +20) 


C3H2 3(2, 2) - 


2(1, 


1) 


< 0.11 








(-30, +30) 


C3H2 3(2, 2) - 


2(1, 


1) 


< 0.15 








(20, -20) 


C3H2 2(1,2) - 


1(0, 


1) 


0.62 (0.20) 


9.9 


(0.6) 


3.3 (1.1) 


(10, -10) 


C3H2 2(1,2) - 


1(0, 


1) 


< 0.11 








(0, 0) 


C3H2 2(1,2) - 


1(0, 


1) 


0.29 (0.04) 


9.3 


(0.1) 


1.1 (0.2) 


(-10, +10) 


C3H2 2(1,2) - 


1(0, 


1) 


0.41 (0.09) 


9.7 


(0.2) 


1.6 (0.6) 


(-20, +20) 


C3H2 2(1,2) - 


1(0, 


1) 


0.25 (0.09) 


9.8 


(0.3) 


1.8 (0.7) 


(-30, +30) 


C3H2 2(1,2) - 


1(0, 


1) 


1.09 (0.17) 


8.4 


(0.4) 


4.6 (0.7) 


(20, -20) 


C2H 1 ^ 






4.77 (0.18) 


10.5 


(0.1) 


3.9 (0.2) 


(10, -10) 


C2H 1 ^ 






3.49 (0.03) 


10.7 


(0.2) 


4.7 (0.3) 
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Table 2 — Continued 



Offset 


Line 




Area 


Vlsr 


AVi/2 


arcsec 






K km 


km 


km s"-*^ 


(0, 0) 


C2H 1 - 


> 


2.32 (0.12) 


9.3 (0.1) 


3.8 (0.3) 


(-10, +10) 


C2H 1^ 


^ 


3.67 (0.10) 


9.3 (0.1) 


5.3 (0.2) 


(-20, +20) 


C2H 1^ 


^ 


3.57 (0.18) 


9.7 (0.1) 


4.5 (0.3) 


(-30, +30) 


C2H 1- 


> 


1.97 (0.23) 


10.3 (0.1) 


2.7 (0.5) 



Note. — Number in parenthesis are 3a errors in each pa- 
rameter. 

^Parameters correspond to the transitions 6(1, 6) — >■ 5(0, 5) 
and 6(0, 6) ^ 5(1, 5) together, which are blended at the same 
frequency. 
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Table 3. LVG results for C3H2 and C2H 



Component 1 Component 2 



Offset 
("/') 


n(H2) 
cm~^ 


iV(C3H2) 

10^2 cm-2 


n(H2) 
cm""^ 


Ar(C3H2) 
10^2 cm-2 


iV(C2H) 
10^^ cm-2 


[C3H2] a 
[C2H] 


(+20, -20) 


<2.4 10^ 


<2.2 


<4.6 10^ 


<l-8 


4.2 


0.010 


(+10, -10) 


8.4 10^ 


1.1 


>3.5 10^ 


4.0 


3.0 


0.017 


(0, 0) 


4.5 10^ 


1.8 


>4.5 10^ 


2.3 


2.0 


0.021 


(-10, +10) 


4.6 10^ 


2.2 


>4.7 10^ 


1.4 


3.2 


0.011 


(-20, +20) 


6.4 10^ 


1.6 


1.1 10^ 


2.4 


3.1 


0.013 


(-30, +30) 


4.5 10^ 


3.1 


4.4 10^ 


1.8 


1.7 


0.029 



Note. — Assumed a kinetic temperature of 50 K. C2H was not separated by com- 
ponents. 

^The abundance ratio involves the whole emission from the two molecules. 
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Table 4. Column densities and relative abundances 



JMoloculo/ riritio 


unit 




Cl(lO.-lO) C2(10.-10) 






C3H2 


10^2 cm- 


-2 


1.1 


4.0 


1.8 


2.3 


C2H 


10^^ cm- 


-2 


1.4 


1.5 


1.2 


0.8 


C0+ 


10^° cm- 


-2 


< 4.8 


< 5.2 


< 4.8 


31.1 


HOC+ 


10^° cm- 


-2 


< 2.6 


9.9 


8.8 


9.6 


HCO 


10^2 cm- 


-2 


< 3.6 


< 3.9 


3.7 


< 1.2 




10^^ cm" 


-2 


9.3 


9.5 


5.9 


5.8 


HC3N 


10^2 cm- 


-2 


2.0 


5.2 


1.1 


4.0 


SiO 


10^° cm" 


-2 


34 


< 22 


11 


< 1.6 




10^° cm" 


-2 


10.9 


14.2 


11.3 


6.6 


HCO+/HOC+ 






> 2700 


930 


830 


450 


CO+/HCO+ 


10-^ 




< 7 


< 6 


< 7 


72 


C3H2/C2H 


10-3 




8 


26 


15 


30 


C3H2/HC3N 






0.6 


0.8 


1.6 


0.6 


HC0/H13C0+ 






< 3.9 


< 4.1 


6.3 


< 2.1 


SiO/Hi3CO+ 






0.37 


< 0.23 


0.19 


< 0.03 



Table 5. Column densities and abundances at selected positions 



Position 
("/') 


A^(HCO) 
10^2 cm-2 


10^2 cm-2 


iV(HC3N) 
10^2 cm-2 


iV(SiO) 
10^1 cin-2 


[HCO] 

TTT I '< 1^1^ -1-1 


[HCO] 
[HC3N] 


[HCO] 
[SiO] 


Comment 


(+12, +24) 


15.8 


1.7 


6.4 


< 2.2 


9.3 


2.5 


> 720 


HCO peak 


(0, +24) 


14.1 


1.0 


4.1 


< 2.2 


14.7 


3.4 


> 640 


HCO peak 


(0, 0) 


< 2.5 


1.6 


11.9 


0.8 


< 1.6 


< 0.2 


< 310 


center 


(+12, -12) 


< 3.7 


1.7 


11.1 


4.3 


< 2.2 


< 0.3 


< 90 


SiO peak 


(+36, +12) 


4.0 


2.8 


18.9 


< 2.0 


1.4 


0.2 


> 200 


Dense envelope 


(-12, -36) 


9.6 


1.8 


25.6 


< 1.6 


5.3 


0.4 


> 600 


Dense envelope 


(0, -12) 


< 3.1 


1.3 


12.0 


6.0 


< 2.4 


< 0.3 


< 50 


SiO peak 



